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Abstract 

We study the r]' — > r/vrvr decay within an effective chiral Lagrangian 
approach in which the lowest lying scalar meson candidates (t(560) and k(900) 
gji I together with the /o(980) and ao(980) are combined into a possible nonet. 

We show that there exists a unique choice of the free parameters of this 
model which, in addition to fitting the vrvr and ttK scattering amplitudes, well 
describes the experimental measurements for the partial decay width of ij' — > 
r/vrvr and the energy dependence of this decay. As a by-product, we estimate 
the ao(980) width to be 70 MeV, in agreement with a new experimental 
analysis. 
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I. INTRODUCTION 



Understanding the status, in general, and the quark content, in particular, of the lowest 
lying scalar mesons is an issue of great current interest. In the cases of the a and the k 
mesons, even their existence has been the subject of many different investigations. One may 
consult refs. [0- [|T6| for a variety of different recent works. 

In the approach upon which this paper is based, a need for a a with a mass around 



560 MeV was found in the analysis of vrvr scattering [T^JThI and a need for a k with a 
mass around 900 MeV was required in order to describe the experimental data on the ttK 



scattering amplitude |T9[. These investigations were carried out in an effective Lagrangian 
framework motivated by the l/N^. approximation to QCD. In this approach, one incorporates 
the contribution of tree Feynman diagrams, computed from a chiral Lagrangian, including 
all possible intermediate states within the energy region of interest. Furthermore, cross- 
ing symmetry is automatic, while the unknown parameters characterizing the scalars are 
adjusted to satisfy the unitarity bounds. Approximate amplitudes satisfying both crossing 
and unitarity are then obtained. For the case of tiK scattering in the I = \ channel the 



analysis of ref. |jT9| may be seen to be consistent with the experimental work of ref. |^0|. The 
experimental analysis characterizes the data by an effective range approximation below 1 



GeV; in the treatment of |T^ it is resolved into the sum of a "current-algebra" piece, vector 
meson exchange pieces and scalar meson exchange pieces. In particular, the presence of a 
K-meson is needed to ensure unitarity. 

Motivated by the evidence for a a and a k, and taking into account other experimentally 
well-established scalars - the /o(980) and the ao(980) - a possible classification of these 
scalars (all below 1 GeV) into a nonet. 
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was studied in ^^ . Since the properties of this scalar nonet are expected to be less stan- 
dard than those of a conventional nonet (like the vectors), the mass piece of the effective 
Lagrangian is allowed to contain extra terms: 

Cmass = -aTT{NN) - bTT{NNM) - cTr(A^)Tr(A^) - dTT{N)TT{NM), (1.2) 

where Ai is the usual quark mass spurion. Retaining just the a and b terms yields "ideal 
mixing" 12^ . The physical particles a and /o which diagonalize the mass matrix are related 



to the basis states N^ and (A^^^ + A'D/v^ by 





COS0, -sin^s , , -^ , /.ON 

smUg cost/s 

where ^^ is the scalar mixing angle. The coefficients a, b, c, and d are determined in terms of 
rricr, rufQ, mao and rrii^, and for a given input set of these masses there are two scalar mixing 
angles. Typical values of the input masses {nia- = 550 MeV, rufo = 980 MeV, niao = 983.5 
MeV and rrif^ = 897 MeV) yield the two possibilities: 

(a) 9s ^ -21°, 

(6) 9s ^ -89°. (1.4) 

In order to determine which of these two possibilities is the correct one, it is necessary to 
study the pattern of scalar-pseudoscalar-pseudoscalar interactions, which are correlated with 
each other by the proposed nonet structure. In this picture, the general form of the SU(3) 
flavor invariant scalar-pseudoscalar-pseudoscalar interaction is: 

Cn4,<^ = Ae'''''edefNX<l)td^4 + BTi (N) Tr {d^<f)d^<f)) + CTr (Nd^cj)) Tr (9^0) 

+ DTr(iV)Tr(9^0)Tr(9^0), (1.5) 

where 0^ [x) is the matrix of the pseudoscalar nonet fields, and A, B, C, D are real parame- 
ters. Derivative coupling to the two pseudoscalars is used to ensure that Eq. (|1.5| ) represents 
the leading term of a chiral invariant expression (see Appendix B of |21|]). It is easy to see 
that all the coupling constants relevant for the study of vrvr and ttK scattering depend only on 
the parameters A and B. The analysis of |2l| then shows that possibility (a) in ( p..4|) for the 



scalar mixing angle is selected as the correct one in the present scheme. The parameters C 
and D were left undetermined in the analysis of [^], as no scalar-pseudoscalar-pseudoscalar 



coupling involving an 7] or rj' was present in the vrvr and nK scattering discussed there. 

In this work we explore the parameter space of C and D in detail by studying the 
rj' -^ rjinT decay, for which there are relatively recent and precise experimental measurements. 
As we will see, the scalar couplings to rj and t]' play a dominant role in the amplitude for 
this decay. 

All the discussion in the present paper will use the same methods and parameters as in 



the previous vrvr and nK scattering papers p8ip!9| . Thus, this work can be thought of as a 



check of that method as well as a test of the basic assumption that the low lying scalars 
are related to each other by belonging to a (broken) flavor SU(3) nonet. In the sense that 
the effective Lagrangian method makes no explicit reference to the quark structure of these 
scalars, the present work may be considered model independent. Note also that only the 



SU(3) flavor structure of tlie scalars is required to construct non-linear cliiral Lagrangians 
describing these interactions |^ . 



"Microscopic" models of low lying scalars have been suggested in which they are variously 



qqqq states in the MIT bag [£J], meson- meson molecules [^ or unitarity corrections due to 
strong meson meson interactions [|l],|12| . All these models involve four quarks and so may be 
related to each other. A "model-independent" effective Lagrangian might be an appropriate 
vehicle for summarizing the common feature of different microscopic models. 

The process rj' —>■ rjTnT has been studied by many authors in chiral symmetric frameworks 
since the early days of "Current-Algebra" . Treatments have used exclusively contact terms 



^ or contact terms plus scalar meson exchanges p0|-|33[|. Ordinarily in the chiral 
perturbation theory approach [^] all effects of resonance exchanges are assumed to be 
"integrated out" and summarized in the complete set of contact terms. However, in the 
case of the 77' (958) decay, the masses of the intermediate a, /o and oq resonances are either 
less than or comparable to 958 MeV. Thus, kinematical dependences due to the propagators 
could be important. The new features of the present treatment include the use of Eq.( |1.2D 
to describe the scalar mesons and mixing angle, the use of Eq.( |1.5| ) to describe the scalar 
coupling constants and a procedure uniform with the discussion of tttt and ttK scattering in 
||T7HT9|. Furthermore, comparison is being made with more recent data. 

This paper is organized as follows. Section II gives our theoretical prediction of the 
7]' -^ rjTnT process as well as the experimental parameterization. The fit to the experiment, 
taking into account the experimental uncertainties, is treated in detail in Section III. Finally, 
Section IV continue a brief summary and discussion. 

II. Tj' -^ r]Tnr DECAY 

Here, we will predict the amplitude for this process in the present model and display the 
experimental data to which it will be compared. 

We assume exact iso-spin invariance which seems consistent with the present experi- 
mental accuracy. The four momenta of the particles are labeled according to the scheme 

V'{p) ~^ vW + ^i(Q'i) + '^2{q2), wherein (7ri,7r2) can stand for either (7r+,7r^) or (7r°,7r°). 
The partial widths are related to the invariant matrix element M{p ^ k + qi + q2) by 

r{r]' -^ r/TT+TT-) = 21(77' -^ ^ttV) = ^— / |M pd$, (2.1) 

2m^/ J 

where the phase space volume element (i$ is 
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FIG. 1. The boundary of integration in Eq. (|2.3| ). 



d^ = {2T,f 6\p -k-q,-q2] 



dk 



dqi 



dc\2 
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(2.2) 



with uo 



mt 



^ + k^ and cUj = J ml. + qi^. After performing the usual phase space integra- 



tion we have 



J- ri'—>r]mT 



1 



GAn^mi 



duJiduj2 I Ml 



The boundary of integration in the ujiuj2 plane for our choice of m„ 
MeV and m^/ = 958 MeV ||3^ is shown in Fig.|l]. 



(2.3) 



137 MeV, m„ = 547 



In the treatment of tttt [|17| , p!8| and t[K [|I^ scattering according to the present approach 
it was found that a reasonable approximation up to the 1 GeV energy range consisted of 
including i) the "current algebra" contact term ii) vector meson tree diagrams and iii) light 
scalar (/o(980), a, k) meson tree diagrams. These were all calculated from a chiral Lagrangian 
with the minimum number of derivatives. For rj' -^ rj-Kix there is a big simplification since 
G-parity conservation shows that no vector meson exchanges are possible at tree level. 
Similarly, the derivative part of the contact term vanishes. 

The individual contributions shown in Fig.E| are then 
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(2.4) 



The total decay amplitude M is the sum of these pieces. The current algebra contribution 
Mc.A. is obtained from the "quark mass" term in the effective Lagrangian (proportional to 
Tr {U + f/^)A^ , where U = exp [2^^/^^). Definitions of the various scalar-pseudoscalar- 
pseudoscalar coupling constants which appear in the a, /o(980) and ao(980) exchange dia- 
grams are given in Appendix A. These involve the coefficients A,B,G,D of Eq. ( |1.2|) ;A 
and B were previously found from tttt and ttK scattering while C and D remain to be 
determined here. The scalar masses are taken as mentioned before Eq.( |1.4| ). Even though 
there is no kinematical possibility for any of the intermediate scalars to be on the "mass 
shell" we include "total width" terms in the propagator denominators in order to agree with 



the previous work ||r7|-p!9|]. The /o and oq exchange terms will be essentially taken to be of 



Breit-Wigner form so G'fQ and G'^q are related to the coupling constants. We take G'fQ = 64.6 



/o 



MeV from ^ ( |§ allows 40-100 MeV) and G'^q = 50 - 100 MeV |3§. The exact value 
of G'^Q will be found from our analysis since it depends on the parameter G. Finally we 
take G'^ = 370 MeV |3^; this is related to a pole position rather than a total Breit-Wigner 
width, a prescription which enables the construction of a vrvr amplitude satisfying both the 
unitarity bounds and crossing symmetry. 

The theoretical expressions in Eq. ( |2.1| )-( P7^ ) will be compared with the experimental 
data on partial decay rates and energy dependence of |M^|. The experimental results for 
the rates are listed [^ as: 



-pexp 
r\-pexp 



0.089 ± 0.010 MeV 
0.084 ±0.012 MeV 



(2.5) 



in agreement with iso-spin invariance. Since we are working in the iso-spin invariant limit 
we will average Q these to obtain: 



T^exp 



0.0872 ± 0.008 MeV, 



(2.6) 



* For the average value x + 5x of measurements Xj + 5xi, we use x 
iJ2i'^i)~^^'^ with the weight Wi = l/{6xi)'^. 
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FIG. 2. Tree Feynman diagrams representing the contributions of (a) the 
current algebra, (b) the a and the /o(980), and (c) the ao(980) terms to the 
decay t( -^ rjnn in our model 



with which the theoretical results will be compared. 



For describing the energy dependence, experimentalists use the Dalitz-like variables p6 
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2 + m^/m^ 2 + m^/m^ 
y = 7:/ (cui + cJaj-lH 77 [m^>-m^) (2.7) 
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with Q = nir 



nir. 



2m7r- As uji and 002 vary over the physical region in Fig.|l|, x ranges 



from about -1.4 to 1.4 and y ranges from -1 to about 1.2. One may expand the matrix 
element, up to an irrelevant overall phase, as 



M = A^'^ 



l + (3iy + p2y^ + I2X' 



+ 



(2.8) 



where A is real while /?!, (32 and 72 are complex. The expansion begins with x^ since M (see 
for example Eq. (|2.4|) ) must be invariant on the interchange gi ^^ q2, which implies x ^^ —x. 
It is found |3^ that this form yields an M^ which fits the experimental data when the yx"^, 
y'^, x^ and y'^x'^ terms are negligible: 



\M\^ = A 



I i2 — 9 

|l + ay| + ex 



+ 



(2.9) 



Here a is complex and c is real. For the decay rj' —>■ rjTT^Tr^, the experimental values are BH 



Re a = -0.058 ±0.013 
Im a = 0.00 ±0.13 
c = 0.00 ±0.03, 



and for the decay 1]' —>■ rjn^ 



TX 



Re a 



-0.08 ±0.03. 



(2.10) 



(2.11) 



As explained before, we compare our results with the average of the experimental data for 
charged and neutral pions. This means we should match our results to 

Re a = -0.062 ±0.012 
Im a = 0.00 ±0.13 

c = 0.00 ±0.03. (2.12) 

The parameter A in Eq. (p.9| ) is determined using Eq. ( |2.6|) . 

Altogether, the experimental data are fit with the four real quantities A, Re a, Im a 
and c. On the other hand the theoretical expression in Eq.( p.4|) is completely fixed if we 



specify just the two real constants C and D in Eq. ( |1.2| ), since everything else is already 
specified. Clearly there is no a priori guarantee that we can fit the data using the present 
model. Furthermore, it is necessary for the expansion of Eq.(|2.4|) to also yield negligible 



higher order terms in Eq. (|2.9| ). We will see in the next section that there in fact exists a 



unique choice of C and D which can fit the experimental data. 

III. FIT TO EXPERIMENT 

Our job is to find the parameters C and D so that |Mp computed from Eq.( |2.4|) agree 
with the experimental form given in Eqs. (|2.9| ) , (p?6| ) and (|2.12|) up to the stated uncer- 
tainties. 

As a preliminary we note that restrictions on the allowed values of C may be obtained 
from experimental information on ao(980) — > vr?] decay. This partial width is given by 

2 

r(ao -^ nri) = ^^^ {mlo -ml- mj) (3.1) 

where q is the center of mass momentum of the final state mesons. Now Eq.( [A10| ) of 
Appendix A shows that 'jcnrri depends on the known values of A and 6p as well as the 



unknown value of C. The Review of Particle Properties |35[ lists the total oq width as 50- 
100 MeV and the irri mode as "dominant". It was estimated in the present model (section 
IV of [^) that r(ao — > KK) is only about 5 MeV so we expect G'^q ~ r(ao ^ tttj) + 5 
MeV. We conservatively expect T{aQ -^ ttt]) to lie in the range 25-100 MeV. This restricts 
C to the two intervals [-21, -13] GeV^^ and [2, 10.5] GeV^^ 

For initial orientation we shall neglect the imaginary terms in the denominators of 
Eq.( |2.4| ). We start by numerically f\ scanning the above two intervals of C and searching 
for the acceptable regions in the CD plane that are consistent with the averaged experi- 
mental partial decay width ( p. 61) . The result of this search is shown in Fig. ^ which also 
shows the analogous intervals when the imaginary terms in Eq.( p.4|) are retained. For C in 
the interval [-13, -21] GeV~^ there is a small acceptable region, whereas for C in [2, 10.5] 
GeV~^ there are two acceptable regions in the form of strips along the C axis. In both 
intervals the thickness of these regions is related to the error in the averaged experimental 
partial decay width in ( p.6|) , and therefore, is the main source of our error estimation in the 
final evaluation of C and D. It turns out that it is a reasonable approximation to neglect 



''^In our computation we choose Og = —20.33°. 
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FIG. 3. Regions consistent with the partial decay width of r/' —>■ rjnn and 
ao(980) -^ TTi]. The semi-closed region on the right is obtained by inclusion 
of the decay widths in the propagators of the intermediate scalars. G'^ = 370 
MeV, G'fQ = 64.4 MeV and G'^q = 100 MeV. The other regions correspond to 
neglecting the decay widths. 



the additional uncertainty associated with the stated error in Re a. In order to further 
restrict the acceptable values of C and D, we compare our predicted energy dependence, 
\M{x, y)\ /|M(0, 0)1 , with the experimental result (|2.9|) and ( p.l2| ) taking Im a = for now 
and c as a fitting parameter. We find that only the region around C = 7 with negative D 
has the required property and therefore we are left with the lower strip in Fig.|^. In Fig.^ 
this region is enlarged; also shown is the line representing a set of "least squared" minima 
on which a is fixed. For a given G, the corresponding minimum is obtained by varying D 
and c. The intersection of this line with the previous region yields the desired G and D 
estimates. Note that the fit improves in the direction of increasing G. The values of G and 
D are displayed in the first column of Table |. 

Now let us include the imaginary terms in the denominators of Eq.( |2.4| ). In our com- 
putation we choose G'^ = 370 MeV and G'^q = 64.6 MeV as were obtained in [|l^, and 
the two extreme possibilities G'^q = 50, 100 MeV. We rescan the GD plane for regions that 
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C5 




FIG. 4. Extracting C and D from two different experimental measurements on 
rj' decay. Circles represent the region consistent with the partial decay width of 
7]' , and the solid line represents the least squared fits of the normalized magnitude 
squared of the decay matrix element to the form (1 + ayY + cx^ with a = 
-0.0615. 



10 



A{GeV-^) 


2.51 


2.51 


2.87 


B{GeV-') 


-1.95 


-1.95 


-2.34 


C{GeV-^) 


7.29 ±0.08 


7.16 ±0.13 


7.25 ±0.10 


D{GeV-^) 


-1.70 ±0.08 


-2.26 ±0.13 


-2.09 ±0.10 


Im a 





-0.12 ±0.27 


-0.16 ±0.20 


c 


-0.004 ±0.031 


-0.014 ±0.033 


-0.013 ±0.033 


S 


1.49 


0.0032 


0.0045 



TABLE I. Extracted parameters from a fit of the normalized magnitude of the 
T)' decay matrix element to the form (1 ± ay)"^ + cx^, with Re a = —0.0615. In 
the first and second columns m,^ = 897 MeV while in the last column m^ = 875 
MeV. The imaginary terms in the propagator denominators were not included 
for column 1. E is the least square deviation with 1701 data points measuring 
the goodness of fit. 



are consistent with the partial decay width ( p.6|) . The result is shown in Fig.^ and is also 
compared with the previous case where no widths were included. This figure shows that in 
the new case, there is no available region for C in the interval [-21, -13] GeV~^. For C in the 
interval [2, 10.5] GeV^^, we have shown in Fig.|^ that the main effect of the inclusion of the 
decay widths is driven by the a width. In Fig.^ we have shown that the uncertainty in G^q 
does not make a substantial difference, in particular in the physical region where C ~ 7. 

We proceed as before, further restricting the available regions in the CD plane by fitting 
the normalized magnitude of the decay matrix element M to the form (|2.9| ) with complex 
a. We set Rea = —0.0615 and fit for Ima and c in this region. We find that the acceptable 
region in this case is very close to the previous region in Fig.|. The result is shown in 
Fig.|^. The two lines correspond to two values of G'^q, and their intersections with the 
acceptable region for t]' partial decay width provide our best points in this plane. We 
however notice that C and D for the value of G'^q = 50 MeV correspond to a value of 
r(ao(980) -^ TCT]) ^ 64 MeV which is greater than the total decay width itself and cannot be 
correct. This consistency check within our computation further restricts the experimentally 
unknown value of G'^q. On the other hand, the intersection of the line corresponding to 



>7r?y 



65 



G'^Q = 100 MeV with the acceptable region of t]' partial decay width gives rao(980)- 

MeV. Therefore we conclude that our computation provides a stable estimate of the partial 

decay width of a0(980) -^ ttt] to be approximately 65 MeV. 

The only other hadronic ag decay mode which has been observed I^Sj is KK; using 
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FIG. 5. The effect of including tlie widtlis in tlie propagators is dominated by 
G'^. In the two parallel regions in the middle, G'^ is removed from its propagator. 
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FIG. 6. The available region consistent with the partial decay width of Vj -^ 
rj-KiT is not sensitive to G"^q in the physical region of C ~ 7. The outer/inner 
regions are obtained with G"^q= 100/50 MeV. 
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FIG. 7. Extracting C and D from two different experimental measurements on 
T)' decay. Circles represent the region which is consistent with the partial decay 
width of 1]' , and lines represent the least squared fits of the normalized magnitude 
of decay matrix element to the form |l-|-a?/| -|- cx^ with Re a = —0.0615, 
G; = 370 MeV and G'jq = 64.6 MeV. 
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r(ao -^ KK) K, 5 MeV |2T| we get an estimate rtoi(ao) ~ 70 MeV. The extracted values of 
C and D and other fitting parameters are hsted in the second column of Table |. Note that 
the goodness of fit improves appreciably when we allow for non-zero widths. 

It is perhaps interesting to display the x and y dependences of our normalized matrix 
element squared |Mp = |M(x,|/)p/|M(0, 0)p. In Fig.^] we show the projections of this two 
dimensional surface onto the y — |Mp and x — |Mp planes. It is clear from the y — |Mp 
projection that |Mp has very little dependence on x. 

The value of the scalar mixing angle dg ~ —21° affects the entire calculation by its 
presence in the formulas [ (|A3|) -( |A17|) ] relating the scalar coupling constants to the parameters 
A^ B^ C and D. Now Og is itself determined by diagonalizing the isoscalar mass squared 
matrix obtained from Eq.( |1.2| ). In this way, 9s depends on the input value of rrif^. The 



value 9s ~ —20.33° corresponds to m^ = 897 MeV but it was shown in |^T| that a range 
865MeV < m^ < 900 MeV gave an acceptable description of ttK scattering. Furthermore, 
reducing m^ to 800 MeV results in the "ideal" case where 9s = 0°. In order to judge the 
sensitivity of our results to changing m^ we repeat the present computation for two lower 
values ttIk = 875 and 800 MeV. As before we scan the CD plane for the acceptable regions 
consistent with the t]' -^ rjinT decay width (|2.6| ). We display the results in Fig.^ which shows 
that the main effect of lowering m^ is in the D > —1 GeV~^ region , far from the physical 
region in which we extract C and D. We see in the same figure that for C ~ 4 — > 8 the effect 
of changing m,^ is negligible. In Fig.lIU] we have displayed these regions together with the 



corresponding least squared fits of the normalized magnitude of the decay matrix element of 
the form (|2.9| ). As we see clearly in this figure, the value of C extracted at the intersection 
of the lines with the strips changes by a very small amount as we go from m«;=897 to 875 
MeV. On the other hand, when we go to the lower value of 772^=800 MeV, the the goodness 
of fit decreases and in particular for C < 7GeV^^ we get unacceptable fits. Furthermore, 
for rrii^ = 800 MeV we get the partial decay width of ao(980) — * tci] to be 124 MeV which is 
greater than the total decay width and is inconsistent. This agrees with the observation in 



| 21| that the values m^ < 875 MeV are not favored. For the value m,^ =875 MeV the details 



of the fit are given in the third column of Table |. 

IV. SUMMARY AND DISCUSSION 

In this work, we studied in detail the rj' —>■ 7/277 decay mode within the framework 
of a model in which the scalar meson candidates cr(560) (discussed in |18[) and /s:(900) 
(discussed in [O) are combined into a nonet together with the /o(980) and the ao(980). 
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squared fits of the normalized magnitudes of decay matrix element to the form 
|1 + ayf + cx^ with Rea = -0.0615. G'^ = 370 MeV and G'jq = 64.6 MeV. 



The scalar mixing angle was calculated ||2T| in terms of these masses using Eq. (|1.2|) and the 



various scalar-pseudoscalar-pseudoscalar coupling constants were calculated in terms of the 
parameters A, B, C and D in Eq.( |1.5D . In the analysis of ref. |]2T]] the parameters A, B 
and 6s were found, but parameters C and D were left undetermined. As rj' decay probes 
these parameters, we have numerically searched this parameter space and found a unique 
C and D which describes the experimental measurements on the partial decay width of the 
rj' -^ rj-Kn as well as its energy dependence. 

Taking into account both the uncertainties in the scalar mixing angle 6s (as reflected in 
the value of m^) and in the rj' — > rjTriT decay width we get for the scalar coupling parameters 

A = 2.51 -^ 2.87 GeV~^ 

B = -1.95 -^ -2.34 GeV"^ 

C = 7.03 -^ 7.39 GeV-^ 

D = -2.39 -^ -1.95 GeV"^ (4.1) 

These numbers are based on combining the second and third columns of Table |I[ The 
coupling constants relevant here are listed in Table |T[ 

As a by-product of the present calculation we obtain an estimate of the ao(980) width 

r (ao(980)) ^ 70MeV (4.2) 

as discussed in Section III. After this work was completed we found a very new experimental 
analysis [^ of the n^p -^ ?/7r+7r^ and n^p -^ rjTT^n reactions which yields the same result 
we have obtained from analysis of the rj' -^ rjinT decay. 

It seems useful to "dissect" our model in order to get a qualitative understanding of 



the rj' -^ rjmr process. Thus we have plotted, in Fig.|TT|, the real and imaginary parts of 
the individual contributions of the terms in Eq.( |2.4| ) to the total decay matrix element. 
These figures again represent projections of the Re M{x, y) and Im M{x, y) surfaces onto 
the lie M — y and Im M — y planes; the small x dependences are thus visible as thickening 
of the curves. First, we observe that the "current-algebra" part of the amplitude, which 
corresponds to the use of the minimal non-linear chiral Lagrangian of pseudoscalar fields, 
is an order of magnitude too small to explain the experimental result by itself. On the 
other hand, the ao(980) exchange contribution is clearly the main one for explaining the 
dominant real part of the amplitude. Nevertheless the other contributions are not negligible. 
For example the cross term 2 [Re M{a)] [Re M(ao)] is of the same order as [Re M(ao)]^. 
Furthermore the a meson exchange is seen to give the largest contribution to Im M for most 
of the kinematical range. 
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Note that we have used just the two input numbers, C and D (over and above the ones 
previously found) to satisfactorily fit the rate and energy distribution of rj' -^ tj-kti. Thus 
in the same framework, with the same parameters, we are explaining tttt |T3| and tiK |T^ 



scattering up to the 1 GeV range as well as rj' —>■ tj-k-k. Our results may then be regarded as 
support for the correctness of both the large Nc approximation motivated approach to low 
energy dynamics being employed as well as the effective Lagrangian model [^ for the low 
lying scalar nonet outlined in the Introduction. Of course, the "microscopic" structure of low 
lying scalars is an interesting puzzle of present day particle physics which seems to require 
a great deal of further experimental and theoretical work for its clarification. For example, 
the study of radiative decays of the 0(1020) is expected [^] to yield useful information. As 



discussed in more detail in |^, the value of the mixing angle 6'^, about —21° and the mass 
spectrum used here are what one would expect with a somewhat distorted form of the qqqq 



model 24 1 . A priori, however, our effective Lagrangian approach can accommodate any 



microscopic model which yields a flavor nonet. 
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APPENDIX A: 

Here we give, for convenience, the explicit form of the scalar-pseudoscalar-pseudoscalar 
interaction [^ . Using isotopic spin invariance, the trilinear Nipcj) interaction from Eq. ( |1.5[ ) 



must have the form: 

- Cn^cI> = -^ (d^Kr ■ df^TTK + h.c?j + -^adf^n ■ d^-rr 

+ "'°^^ df,KT ■ aod^K + 'j^kv {^.0^X8^,7] + h.c.) + '^^Kr,' {Rd^Kd^r]' + h.c.) 
+ 7ao7rr?ao " d^Tzd^rj + 7ao^v^o ■ d^nd^T]' 

+ lfomfo9,,vd^V + lfovv'fo9,,vd^v' + Ifov'v'fo^t^v'dt^v', (Al) 

where the 7's are the coupling constants. The fields which appear in this expression are the 
isomultiplets: 
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K 



K+ 






K 



, H = (k K°) 



n 



V2 



l7ri=F«7r2J, 71 =7l3, 



V2 



(aoi =F iao2) 



% — C^OS) 



(A2) 



in addition to the isosinglets a, /o, 77 and rj'. The 7's are related to parameters A, B, C, D 
of Eq. dlj) by 

(A3) 
(A4) 
(A5) 
(A6) 
(A7) 
(A8) 
(A9) 
(AlO) 
(All) 



IkKit = laoKK = -2v4 

'janw = 2Bsm9s - V2{B - A)cos9s 
j^KK = 2(25 - v4)sin^, - 2V2Bcos9, 
7^^^^ = ^/2{A - B)smes - 2ficos6's 
jj^j^j^ = 2{A - 25)cos^, - 2V25sin^, 
IkKt) = CsinOp — v2(C — A)cos6p 
'IkKti' = y2{A — C)sin6'p — CcosOp 
lao-KTi = (C - 2y4)sin6'p - \/2Ccos6p 
lao-KT)' = (2 A — C)cos9p — v2Csin^p 



1 



7. 



crrj'rj' 



V2{B + D)--{C + 2A + 4D)sm29p + V2{C + D)cos^ep 

{B + D)--^{C + 2D)sm29p + {A + D)cos^9p + Csin^Op 
V2{B + D) + -{C + 2A + 4:D)sm2ep + V2{C + D)sm'^9p 

{B + D) + ^{C + 2D)sm26p + {A + D)sm'^ep + Ccos^^p 

v2 



sin6's 

cos6's 
sin^s 

0086*, 



7. 



arjri' 



[V2{C + D)sm29p + {C + 2A + AD)cos2ep 
- [V2{C + 2D)cos20p +{A-C + D)sm2e. 



sin6's 

0086*, 



Ifovv 



Ifov'v' 



V2{B + D) + -{C + 2A + AD)sm2ep - V2{C + D)cos^ep 

{B + D)-^{C + 2D)sm2ep + {A + D)cos% + Csin^^p 

v2 

V2{B + D) + -{C + 2A + AD)sm2ep + V2{C + D)sin^9p 

{B + D) + ^{C + 2D)sm29p + {A + D)sm'^9p + Ccos^dp 

V2 



cos6's 
sin^s 

COS^s 

sin6'o 



7/, 



ow 



- [V2{C + L))sin20p + (C + 2A + AD)co?,2dp 
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COS^s 



(A12) 



(A13) 



(A14) 



(A15) 



(A16) 



Tttttt 


7.27 


7.27 


8.36 


laKK 


9.63 


9.63 


10.44 


la-qri 


3.90 


4.11 


4.30 


^ariT)' 


1.25 


2.65 


2.61 


''larfr)' 


-3.82 


-1.43 


-2.09 


T/TTTT 


1.47 


1.47 


2.53 


ISKK 


10.11 


10.11 


12.76 


^fm 


1.50 


1.72 


2.78 


Ifm' 


-10.19 


-9.01 


-9.34 


Ifv'v' 


1.04 


2.60 


2.04 


7a7rT) 


-6.87 


-6.80 


-7.28 


^anr]' 


-8.02 


-7.80 


-7.38 



TABLE II. Predicted coupling constants corresponding to the columns in Table 
ID. All units are in GeV^^. 



\/2(C + 2D)cos2ep +{A-C + D)sm29^ 



sin0„ 



(A17) 



where 6s is the scalar mixing angle defined in Eq. (|1.3|) while 6p is the pseudoscalar mixing 
angle defined by 



V 
V' 



cosOp —sinOp 
sinOp cos9p 



l)/V2' 



(A18) 



where rj and rj' are the fields which diagonalize the pseudoscalar squared mass matrix. We 
adopt here the conventional value 9p ^ 37°. (see |21] for additional discussion.) 
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